Complex and real unconventional Bose-Einstein condensations in high orbital bands 
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We perform the theoretical study on the unconventional Bose-Einstein condensations (UBEC) in 
the high bands of optical lattices observed by Hemmerich's group. These exotic states are char- 
acterized by complex-valued condensate wavefunctions with nodal points, or real-valued ones with 
nodal lines, thus are beyond the "no-node" paradigm of the conventional BECs. A quantum phase 
transition is driven by the competition between the single particle band and interaction energies. 
The complex UBECs spontaneously break time-reversal symmetry, exhibiting a vortex-antivortex 
lattice structure. 



PACS numbers: 03.75.Nt, 03.75.Lm, 05.30.Jp, 05.30.Rt 



Quantum wavefunctions are generally complex- valued. 
However, the usual ground state wavefunctions of bosons 
are very restricted because they are positive-definite as 
stated in the "no-node" theorem This theorem ap- 
plies under very general conditions: the kinetic energy 
is unfrustrated (e.g. the Laplacian-type) ; the single par- 
ticle potential can be arbitrary; the two-body interac- 
tion depends only on coordinates. Mathematically, it 
is a direct consequence of the Perron- Frobenius theorem 
of matrix analysis |2j. This theorem implies that time- 
reversal (TR) symmetry cannot be spontaneously broken 
in various ground states of bosons, including supcrfluid, 
Mott-insulating, and supersolid states. 

The "no-node" theorem, however, is a ground state 
property, thus it does not apply to meta-stable excited 
states of bosons. This opens up a possibility for "un- 
conventional" states of bosons beyond the "no- node" 
paradigm Similarly to unconventional supercon- 

ductors, in unconventional Bose-Einstein condensations 
(UBEC), the condensate wavefunctions form non-trivial 
representations of the lattice symmetry groups. However, 
a major difference exists. Cooper pairs have the center of 
mass motion and the relative motion between two elec- 
trons of the pair. In unconventional superconductors, it 
is the relative motion that is non-trivial. The degree of 
freedom of the relative motion does not exist in the sin- 
gle boson BEC. In UBECs, the condensate wavefunctions 
are non-trivial. 

Considerable efforts have been made to study uncon- 
ventional states of bosons both experimentally and theo- 
retically. Among the most exciting achievements are the 
realizations of the meta-stable excited states of bosons 
in high orbital bands 0-0] > which leads to the opportu- 
nity to the study of the UBECs @-[I3, and other exotic 
properties |16h21| . Below are some recent experimen- 
tal results. Sebby-Strabley et al. succeeded in pump- 
ing a large fraction of bosons into the excited bands in 
a double- well lattice [J]. Mueller et al. observed the 
quasi-lD phase coherence pattern by exciting bosons into 
the p-orbital bands in the cubic lattice [5[. An impor- 
tant progress was made by the group of Hcmmerich Q : 
the UBECs in the sp-hybridized orbital bands were real- 
ized in a checkerboard-like lattice, which allows to estab- 
lish the fully cross-dimensional coherence. More recently, 



UBECs in even higher orbital bands have been observed 
in the same group Q- 

In this paper, we present the theoretical study on 
UBECs observed in the second, or, the first excited band, 
of the checkerboard optical lattice. This band is of a hy- 
bridized nature between the s-orbital of the shallower 
sites and the p-orbitals of the deeper sites. The lattice 
asymmetry favors a real- valued condensate wavefunction 
with nodal lines, while interactions favor a complex- 
valued one with nodal points. By solving the Gross- 
Pitaevskii (GP) equation for these meta-stable conden- 
sates, we find that tuning the lattice asymmetry drives 
the phase transition between these two types of UBECs 
in a good agreement with experimental observations. 

We introduce the optical lattice employed in the ex- 
periment @. Each unit cell consists of two sites with 
different depths (denoted A and B below) as shown in 
Fig. Q] (a). (A similar lattice potential with different 
parameters has been plotted in Ref. &]). The lattice is 
constructed by the interference pattern of phase coher- 
ent laser beams along ±x and indirections generated 
from a single laser through beam splitters and reflectors. 
The optical potential reads 
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where y and i are unit vectors describing light polar- 
izations; ki is the laser wavevector; e < 1 and rj < 1 
describe the imperfect reflection and transmission effi- 
ciencies; 9 is the phase difference between beams along x 
and y-directions; a is used to tune the lattice asymmetry 
by rotating the light polarization out of the ^-direction. 

The point group symmetry of this lattice is analyzed 
below. We start from the ideal case of e = 1 with 
a = 0° and 9 = 90°, at which A and £?-sites are equiv- 
alent. At 77 < 1, the lattice has the reflection sym- 
metries with respect to both the x and y-axes, thus 
the lattice is orthorhombic. Next we keep e = 1 and 
a = 0° but set 9 away from 90°. Then the unit 
cell includes both A and B-sites. The primitive lat- 
tice vectors are ao(e x ± e y ) where ao = n/ki as shown 
in Fig. Q] (b). The optical potential becomes V = 
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FIG. 1: (a) The optical lattice with the reflection symmetry 
with respect to the x-axis and the parameter values: r\ — 0.95, 
e = 0.81, 9 = 95.4°, a = a = 36° and Vo = 6.2E r . The A- 
sites have deeper potential depth than those of B-sites. (b) 
The basis vectors of the double-well lattice. 




FIG. 2: (a) The energy spectra for the second band, the pa- 
rameter values are the same as Fig[T] except a — 0°. (b) The 
spectra of (a) along the line from (0, ~) to (f , 0). 



— ■^■(cos2fc;a; + rj 2 cos 2kiy + 4r] cos 8 cos k[X cos kiy) . 8 
controls the potential difference between A and B sites. 
The point group symmetry remains orthorhombic. Now 
we move to the realistic case of e < 1. The unit cell re- 
mains double-well-shaped and the primitive lattice vec- 
tors are the same. However, the orthorhombic symmetry 
is broken and there is no point group symmetry for gen- 
eral values of parameters. This asymmetry can be par- 
tially compensated by setting ao = cos _1 e. We denote 
this configuration as "symmetric" and other ones with 
a ^ ao as "asymmetric" below. The symmetric lattice 
potential becomes V = — -y-e(e cos 2kix + rj 2 cos 2kiy) — 
Vorjecoskix[cos(kiy + 8) + e 2 cos(kiy~8)] , which has the 
reflection symmetry with respect to the £-axis but not to 
the y- axis. 

Next we calculate the band structures. The reciprocal 
lattice vectors are defined as G myn — mb\ + nb 2 with 
&x,2 = (if ! -)• The single particle Hamiltonian reads as 
H = -h 2 V 2 /(2M) + V(r) where M is the boson mass. 
Using the plane wave basis, the diagonal matrix elements 
are (k + G mn \H \k + G mn ) = E r {[ak x /w + (m - n)] 2 + 
[aky/n + (to + n)} 2 }, where E r = H 2 tt 2 / (2Ma 2 ) is the 
recoil energy. The off-diagonal matrix elements read 

(k\H \k + G ±lfi ) = — jmicosae^ 6 +e ±l6 ), 

(k\H \k + G ,±i) = — jr;(cosae ±ie + e 2 e^), 

- Vo 
(k\H \k + G±i, T i) = — — ecosa, 

(k\H \k + G ±1 , ±1 ) = -^eri 2 cosa. (2) 

We focus on the second band into which bosons are 
pumped [6|. There are four points in the Brillouin 
zone (BZ), i.e., O = (0,0), K 1<2 = (i^i), and 
X = (—,—), at which the Bloch wavefunctions are TR 
invariant, and thus real- valued. The band spectra are 
symmetric with respect to these points, which means that 
they are local energy extrema or saddle points. For the 
symmetric lattice with a = ao , the second band has dou- 
bly degenerate energy minima of the states ipKt 2 located 
at K12, respectively. For the asymmetric case, the degen- 
eracy between and ipK 2 are lifted. For a < (>)cvq, 



Ki(K 2 ) become the band minimum, respectively. The 
energy spectra of a — is shown in Fig. [5] (a) and (b) 
(A similar energy spectrum with different parameters has 
been plotted in Ref.[6j). 

The real space distributions of tpKi 2 are a ^ so calcu- 
lated. Their nodal lines pass the centers of the deeper 
sites of A. Thus the orbital component on the ^4-sites is 
of the p-type and that on the shallower sites of B is of the 
s-type. In fact, the p-orbital configurations of 2 (^) m 
the A-sites are actually not exactly along the directions 
of e x ± e y because of the lack of the tetragonal symmetry. 
This point is mostly clear in the case of strong potentials 
so that we can define local orbitals on each site. Even for 
the symmetric lattice, the p x and p y -orbitals on the A- 
sites can be defined according to their parities under the 
reflection with respect to the x-axis. However, they are 
non-degenerate. The orbital components of ipKi 2 (r) are 
nearly the same on A-sites, i.e., mostly the lower energy 
p-orbital slightly hybridized with the higher one. The or- 
thogonality of these two states comes from their different 
lattice momenta. 

Interactions determine the configurations of UBECs in 
the presence of degenerate band minima. Any linear 
superposition among them gives rise to the condensate 
wavefunctions with the same kinetic energy. However, 
interactions break this degeneracy. Previous studies on 
p-orbital BECs based on tight-binding models predicted 
linear superpositions between two Bloch wavefunctions 
at degenerate band minima with a phase difference ±5. 
Such a condensate breaks TR symmetry spontaneously 
d, H . Bosons on p-orbital sites aggregate into the p x ± ip y 
states to reduce their repulsive interaction energy. This 
is a result of the second Hund's rule: complex p-orbitals 
are spatially more extended than the real orbitals, and 
thus bosons have more room to avoid each other. 

The optical potential in the current experiment is shal- 
low, thus the system is in the weak correlation regime [f|. 
Instead of the tight-binding model, we use the GP equa- 
tion. Because of the absence of the lattice potential along 
the z-axis, we neglect the z-dependence of the conden- 
sate wavefunction. We only consider its distribution &(f) 
in the icy-plane. It is normalized as ^ J d 2 r\^{f)\ 2 — 1 
where f d 2 r integrates over one unit cell with the area 
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FIG. 3: The distributions of (a) the phase and (b) density 
patterns of the complex UBEC. The parameter values are the 
same as in Fig[T]except gpo = 0.6i5 r and a = 36°. The vortex 
and anti- vortex cores are located in the centers of A-sites. 

of n = 2oq. The GP equation reads 

{ ~ W + V{f) + 5Po|*W| 2 }*(r) = (3) 

where po = Nq/V is the average 3D density with No the 
total boson number in the condensate and V is the 3D 
volume of the system; g is the s-wave scattering interac- 
tion parameter. In the calculations below, various values 
of interaction parameters gpo are used from up to E r . 

(a) (b) (c) 




FIG. 4: Density distribution in the time-of-flight spectrum 
for (a)complex condensate in the symmetric case (a = 36.0°) 

(b) complex condensate in the asymmetric case (a = 35.5°). 

(c) real condensate (a = 34.5°). Other parameters values are 
the same as Fig[T] except a. 

Although Eq. [3J looks the same as the usual GP equa- 
tion, the marked difference is that ^(r) is not the ground 
state condensate but the meta-stable one belonging to 
the second band. The non-linearity of the GP equation 
allows mixing between different Bloch wave states. Let 
us start from the symmetric lattice with a = otQ. Eq. 
[3J is solved self-consistently as follows. We define the 
renormalized potential as V e ff(f) = V(f) + gpo|^E'(r)| 2 , 
and solve the corresponding renormalized band structure. 
Then the condensate wavefunction is optimized to min- 
imize the total energy, which in turn determines V e //. 
The renormalized band structure is similar to the free 
one, which still has two degenerate band minima at K\ 2- 
We define the condensate wavefunction as 

^(r) = cos 8 ip Kl (r) + e i0 sin 8 ip Kz (r). (4) 

The total energy reaches minimal at 8 = ? and cf> — 
±7^. These complex condensate wavefunctions only have 
nodal points, while the real ones ipK t 2 have nodal lines. 
The complex ones are spatially more uniform, and thus 



are favored by interactions. We plot the phase and den- 
sity patterns of this condensate in Fig. [3J which exhibit 
a vortex-antivortex lattice structure. The vortex and 
anti-vortex cores are located alternatively at centers of 
A-sites, at which the antiferromagnetic order of orbital 
angular momentum develops. For every closest four B- 
sites, their phases wind around the central A-site follow- 
ing the same vorticity. This is similar to the case of the 
tight-binding models 3, 8]. The Bragg peaks in the time 
of flight (TOF) spectra are located at (m + 3)61 + nb 2 
and mbi + (n + -5)^2 as observed in the experiment @. In 
particular, the four peaks of ±5^1,2 = (i^o - ' ^So - ) are 
strongest with equal intensities, as shown in Fig|4] (a). 
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FIG. 5: (a)The phase diagram as a and the interaction 
strength gpo- Other parameters values are the same as FigO 
(b) The condensate fractions of ipK x in the complex UBEC 
$ = cos S ipfd ± i sin S i[>k 2 , the parameter values are the 
same as FigQ] except a and gpo = 0.6i5 r . 

Now we move to the asymmetric lattice whose free 
band structure minimum is non-degenerate. The com- 
plex condensates are favored by interactions, and thus 
should be stable at asymmetries weak enough. Certainly, 
at large asymmetries, the real condensate wins due to the 
gain of band energy. This picture is explicitly confirmed 
by the phase diagram calculated by GP equation. As 
shown in Fig. [5] (a), for a given value of the interaction 
strength gpo, the complex condensate in the form of Eq. 
2] is stable in a finite parameter range from a\ to 0:2, be- 
yond this regime the condensate changes to the real one, 
and the TOF spectra of such a real condensate only con- 
tain peaks of (m + \)b\ or (m + ^)? 2 , as shown in FigfJ] 
(c) 

In the complex condensate, the relative phase <f> be- 
tween ipKi, 3 is always ±-|, i.e., ^ and ^* are degenerate 
as TR partners; 8 is asymmetry dependent. The spatial 
asymmetry of ^(r)! 2 depends on that of the bare po- 
tential V. However, V e ff, a combination of V and l^l 2 , 
becomes symmetric. Without loss of generality, \l/(r) is 
expanded in terms of two orthonormal real wavefunctions 
^1,2 (j 5 ) in the same way as in Eq. [4] by replacing ipKi 2 
with "01,2- Apparently, both fy(r) and , F*(r) satisfy Eq. 
[3J and yield the same V e f f . The corresponding renormal- 
ized single particle Hamiltonian, — ft 2 V 2 /(2M) + V e ff, 
has degenerate band minima ^12- However, please note 
that the superposition principle does not apply to the 
non-linear GP equation: ipifi are not solutions to Eq. [3J 
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^(r)! 2 is also asymmetric depending on the asymmetry 
of the bare potential V. The TOF spectra still exhibit 
four dominant peaks at ±(^, ^) and ±(-5^, as 
shown in Fig [J] (b). The relative intensities of these two 
pairs of peaks depend on the lattice asymmetry, which 
can be reflected by the condensation fractions ipKi in the 
complex condensate, as plotted in Fig. [5](b). An obser- 
vation of the asymmetric peaks at ±i6 12 at ct\ < a < a 2 
would provide a supporting evidence for the complex con- 
densates. The TOF spectra lack phase information, thus 
the observation of the symmetric peaks ±561,2 at cyq [6| 
could be interpreted as the phase separation of real con- 
densates of ipK-i 2 , or an incoherent mixing between them. 
However, in these scenarios, the lattice asymmetry lifts 
the degeneracy and only leads to one pair of peaks. Even 
two condensates could coexist forming domains, their na- 
ture is of hysteresis. The condensate fraction of ipK 1 in 
the complex condensate should not follow that plotted in 
Fig. Eh)- 



(a) r, a = %i (b) 11 = 




FIG. 6: Phase diagram as a function of rx,T2 predicted by 
EqUfor (a) g > 0; (b) g < 0. 

For a better understanding of phase transitions 
between real and complex UBECs, we construct a 
Ginzburg-Landau (GL) free energy as: 

F = -ri|* Kl | 2 -r 2 |*K 2 | 2 +5i|*Ki| 4 + 52|*if 2 | 4 
+ g 3 \^ Kl \ 2 \^ K2 \ 2 + gi (^* Kl ^* Kl ^K 2 ^K 2 +h.c), (5) 

where ^K lt2 = ^1,2 e* 91 - 2 describe the condensate order 
parameters at K\ft #1,2 are the phases of the condensates 
of ^k 1<2 an d ^1,2 are real as explained before. Although 
^k 12 do not couple at the quadratic level due to the re- 
quirement of translational symmetry, they do couple at 
the quartic level as in the g\ term because ±2{K\ — K2) 
equals reciprocal lattice vectors. $4 is positive for repul- 
sive interactions which favors the relative phase differ- 



ence #i — 62 = ±§ , thus the free energy in Eq.© can be 
reduced to: 

F = -ril^K, -r 2 ^K 2 +9l^Ki +92^K 2 +9^K 1 i > K 2 , (6) 

in which g = g 3 — 2g^. We define G = ^9\92 — g 2 and 

9it92tG > as required by the thermodynamic stability 

condition. In the superfluid regime, the complex UBEC 

is characterized by the non-zero values of both ^k 12 , 

while the real BECs correspond to one of these values 

being zero. Without loss of generality, we fix gi,g2,g 

and plot the phase diagram of the superfluid regime as 

a function of r\ 1 r2- As shown in Figj6l for g > 0, the 

complex BECs occur when < — < — , where both 

. . . . 92 . T2 9 . 

ri,r2 are positive. It is interesting to notice that for 

g < 0, the complex BEC can exist even one channel is 
off-critical (ri < or r 2 < 0), which means that in this 
case, the complex BEC is purely induced by interaction. 

As interaction increases, and the system is brought into 
the Mott insulating regime. Nevertheless, at least in the 
weakly insulating regime, the suppress of the superfluid- 
ity ordering is due to phase fluctuations, and the magni- 
tudes of l^Kt 2 1 remain nonzero. Though B\ and 9 2 are 
disordered such that (^k 1i2 ) — 0, their relative phase 
81 — 62 = ±§. This indicates a TR breaking order with 
a bilinear form of ^ Kl , 2 as L = i(^* Kl ^K 2 - *k 2 *Ki) 
in the Mott insulating state. Its physical meaning here 
remains the staggered circulating currents, i.e., this ex- 
otic Mott insulating states preserve the antiferromagnetic 
OAM order of the complex BECs but not the global phase 
coherence. 

In summary, we have studied the UBECs observed 
in high orbitals bands in Ref.[6|. The unconventional 
condensate wavefunctions can be real and TR invari- 
ant with nodal lines, or complex breaking TR symmetry 
with nodal points. In both cases, translational symmetry 
is broken due to the nonzero condensation wavevectors, 
thus these UBECs can be considered as unconventional 
supersolid states. The interplay between lattice asymme- 
try and interactions drives the transition between them. 
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